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Abstract 

The d-dimensional string generated gravity models lead to Einstein- 
Maxwell equations with quadratic order correction term called the 
Gauss-Bonnet term. We calculate the quasinormal modes for the d- 
dimensional charged black hole in the framework of this model. The 
quasinormal spectrum essentially depends upon the Gauss-Bonnet cou- 
pling parameter a which is related to the string scale, and is totally 
different from that for black holes derived from Einstein action. In 
particular, at large a the quasinormal modes are proportional to a, 
and thus, have greater oscillation frequency and damping rate, while 
as a goes to zero the qusinormal modes approach their Schwarzschild 
values. In contrary to Einstein theory black hole behavior, the damp- 
ing rate of the charged GB black hole as a function of charge does not 
contain a chracteristic maximum, but instead the monotonic falling 
down is observed. In addition, there have been obtained an asymp- 
totic formula for large multipole numbers. 



1 Introduction 



Black holes in greater than four space-time dimensions have attracted a lot of 
interest within different theories. One of them is the supersymmetric string 
theory: its slope expansion yields corrections to the Einstein action. These 
corrections contain higher then the first powers of curvature. The quadratic 
term in curvature is the leading correcting term and can affect the graviton 
excitation spectrum near flat space. Thus one can consider the low-energy 
limit of string theories as an effective models of gravity in higher dimensions. 
The full low-energy theory must solve the problem of singularities in general 
relativity. The so called Gauss-Bonnet invariant is one of the most promising 
candidates for such quadratic curvature corrections to the Einstein action. 

In addition, Lovelock came to this action when investigating the fol- 
lowing question in classical theory of gravitation: what is the most general 
divergence free second rank symmetric tensor in D-dimensional space-time, 
constructed form the metric and its derivatives up to the second order?. 

The spherically symmetric solution describing neutral static black hole 
in Gauss-Bonnet gravity was obtained by Boulwar and Deser in pp. The 
charged black hole solution was found by Wiltshire [2]. Since that time, 
the properties of Gauss-Bonnet black holes were studied in a lot of papers 
[H]-jl5j. with special interest to black holes in anti-de Sitter space-time as 
that having interpretation due to ADS/CFT correspondence. 

As is known, one of the most important characteristic of a black hole, 
which play a fundamental role for both classical (see ^B] f° r a review) and 
quantum description is its quasinormal spectrum. The quasinormal 
modes are resonance oscillations which dominate at late times of the black 
hole response to an external perturbation. Recently a lot of papers (see for 
instance |18j - |35| and references therein) have been devoted to this subject 
motivated either by observational aspects or quantum aspects coming from 
Loop Quantum Gravity and ADS/CFT correspondence. 

At the same time the higher dimensional black holes have been actively 
investigated recently within different extra dimensional scenarios (see 
for a recent review). There, it is expected that the Planck scale is of order 
of TeV, and thereby the quantum gravity may show itself at LHC where 
it is expected to observe mini-black holes. The Gauss-Bonnet black holes 
are being studied in this context as well (see |H] and references therein). 
In this connection, the quasinormal modes which are independent of the 
way of their excitation, and, thereby represent the "footprints" of a black 
hole, could give us the footprints" of the Gauss-Bonnet coupling, i.e. of low 
energy regime of quantum gravity. 
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Therefore we are interested here, what will happen with quasinormal 
spectrum of the black hole if the Gauss-Bonnet correction term will be taken 
into consideration. To our surprise the first, and as far as we are aware 
the only, paper on this subject was written yet 15 years ago [33]. There, 
some modes for few values of the slope parameter were evaluated when 
considering the scalar wave in the five and sixth dimensional neutral Gauss- 
Bonnet black hole background. Here we are performing a thorough study 
of the Gauss-Bonnet black hole quasinormal modes for both neutral and 
charged black hole which include the whole range of the slope parameter 
a, and different multipole numbers I and different dimensionality of space- 
time D. In particular, we investigate the asymptotic regime of large a 
and I. It is shown that the quasinormal behavior is totally different from 
that for the ordinary black hole corresponding to Einstein action. Thus, 
the quasinormal behavior is crucially depend upon the slope parameter a: 
the QN modes approach their Schwarzschild values when a goes to zero, 
and are proportional to a at large values of a. Thus, influence of Gauss- 
Bonnet (GB) coupling lead to increasing of the oscillation frequencies and 
to a greater decay rate in the regime of large GB coupling a. For moderate 
values of a, the the damping rate, as a function of a, is decreasing reaches 
some minimum, and then increasing as ~ a. 

The dependence of the spectrum on the charge of the black hole is also 
different from that for Reissner-Nordstrom black hole: the imaginary part is 
just monotonically decreasing as the charge is increased, while for charged 
black holes derived from Einstein action such as Reissner-Nordstom black 
hole or its dilaton extension, one has the characteristic maximum of the 
imaginary part of to somewhere near the 0.8 of its extremal value of charge. 

The paper is organized as follows: in Sec. II we give some preliminary 
information on Gauss-Bonnet black hole perturbations. Section III deals 
with calculation of the quasinormal spectrum for Gauss-Bonnet black hole. 
Sec IV is devoted to summary of the result and some future perspectives of 
this work. 

2 Basic formulas 

The Lagrangian for the Einstein action and the Gauss-Bonnet action in the 
D-dimensional space-time model has the form: 
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+ ai j d D xV=g(RabcdR abcd - 4R cd R cd + R 2 ) (1) 

where ai is (positive) Gauss-Bonnet coupling constant which is related to 
the Regge slope parameter or string scale. 

The black hole solution for this action is described by the metric: 

ds 2 = -f{r)dt 2 + r l {r)dr 2 + r 2 dn 2 D _ 2 . (2) 

Here the asymptotically flat solution corresponds to the following form 
of the function / [""]: 



Here \x is the integration constant which is related to the black hole mass, 
and, a and at are connected by the relation: 

a = 1&kGd(D — 3)(.D — 4)a/ (4) 

Note that in four space-time dimensions the Gauss-Bonnet term is the 
topological invariant "3J"|. The gravitational and thermodynamic properties 
of the above black holes were studied in a lot of papers recently [3J, [I], 
[TU] . [TT] . [T""] . "14"] . [1*5] . including the extension to the non-zero lambda 
term case. There are some range of parameter a and \i for which a black 
hole shows either thermodynamical or gravitational ( or both) instability 
jlUj . |14j . Yet for D > 6 and D = 5 the Gauss-Bonnet black hole is stable 
against tensor type of gravitational perturbations "TH] , which, as was shown 
by Gibbons and Hartnoll "37] is, in a sense, the decisive type of gravitational 
perturbations which produce instability at least in the Einstein theory black 
holes. The charged Gauss-Bonnet black hole corresponds to the following 
form of the function / [5] : 



M _ f_ _ f_ \ 8afi 4aQ 2 

} [r} ~ + 2^ " 2^ V r D ~l ~ 2tt(D - 2)(D - 3)r 2D " 4 ' (; " 



Here Q is the charge of the black hole. In the latter case, there will be 
timelike singularity shielded by two horizons which tend to coincide as the 
charge Q goes to its extremal values Q ex t- The extremal value of charge is 
determined by the relation: 

2(D-3) D-5 2 (D-4) _ Ql xt _ q (Q) 

r ext + a D _ 3 r ext 2tt(D-2)(D-3)~ U W 
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where 



2 \ 1/2 



& - - 5), + (\(D - 5)V + ^g:^ 3) j • P) 

To consider scalar perturbations in the background of the above Gauss- 
Bonnet black hole we need to consider the test scalar field in the GB black 
hole background, i.e the consider the Klein-Gordon-Fock equartion: 

□$ = ^L( 5 ^x/=ff*, M )„ = 0. (8) 

V-g 

After the change of the wave function <!>, and going over to the tortoise co- 
ordinate dr* = dr/f(r), the last equation can be reduced to the Schrodinger 
wave-like form: 

W - F(r*)) * = 0. (9) 
An effective potential has the form 

- m ( (p - 2 » r ' 2 p - V ) + + SLtgzj)) po) 

The tortoise coordinate r* is defined on the interval (— oo, +oo) in such a 
way that the spatial infinity r = +oo corresponds to r* = oo, while the event 
horizon corresponds to r* = — oo. The above effective potential is positively 
defined and has the form of the potential barrier which approaches constant 
values at both spatial infinity and event horizon. 



3 Quasinormal modes of GB black hole 

Under the choice of the positive sign of the real part of u (u = Reu> — Imuj), 
QNMs satisfy the following boundary conditions 

*(r*) ~ C±ex.p(±iujr*), r — > ±oo, (11) 

corresponding to purely in-going waves at the event horizon and purely out- 
going waves at infinity. 

To find the quasinormal modes of the black hole whose the effective 
potential has the form of potential barrier like that of the formula (9), one 
can use the WKB approach. The latter was first used for calculations of 
quasinormal modes in the first order by Schutz and Will [HE] , and extended 
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by Iyer and Will to the third order beyond the eikonal approximation |39| . 
and then extended to sixth WKB order in [IUJ. The WKB approach for 
finding QN modes is used in a lot of papers, and in particular the WKB 
corrections up to 6th order have been used recently in [13], [42] . The 
WKB formula has the form [10] : 

* U J r^TTU ~ L 2 ~ L 3 - L 4 - L 5 - L 6 = n + ^, (12) 

where Vq is the height and Vq is the second derivative with respect to the 
tortoise coordinate of the potential at the maximum. L 2 and L3 can be 
found in [39] . L4, L5 and are presented in j^Oj; the corrections depend 
on the value of the potential and higher derivatives of it at the maximum. 
Below we shall discuss the obtained qusinormal modes classified according 
to the multipole number for charged and neutral GB black hole. 

From here and on we shall take [i = 1 for any dimension D, thus, in 
different space-time dimensions the units are different. 



3.1 The neutral black hole, 1=0 

Within the WKB approach, we shall use here, I = modes represent an 
exclusive case, since the unmodified WKB method is a good approximation 
when I > n, being not so good for I = n (and in fact being inapplicable when 
I > n). As is known from comparison with numerical computations jHH] ; the 
error at I = scalar modes can reach even 10 per cent for the Schwarzshild 
black hole case. So one should be careful when interpreting the WKB results 
for this case. The I = modes for cZ = 5,6, 7, 8-dimensional black holes are 
represented in the following tables for different values of a : 
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a [JJ = 5) 


i ~> 

Ke lo 


-1m lo 


a {V = 7) 


Ke lo 


-1m lo 


1/10 


0.389935 


0.256159 


1/10 


1.11738 


0.546056 


1/5 


0.396034 


0.250548 


1/5 


1.13699 


0.529543 


1/2 


0.429741 


0.208293 


1/2 


1.29469 


0.395111 


5 






5 


1.39823 


0.574472 


10 






10 


1.48053 


0.385616 


20 






20 


2.00368 


0.56458 


a (£> = 6) 


Re lo 


-Im lo 


a (D = 8) 


Re lo 


-Im lo 


1/10 


0.735854 


0.402416 


1/10 


1.51702 


0.694245 


1/5 


0.748053 


0.391049 


1/5 


1.54463 


0.673566 


1/2 


0.83530 


0.304837 


1/2 


1.7854 


0.488086 


5 


0.906661 


0.144820 


5 


1.47941 


0.868859 


10 


1.513624 


0.456935 


10 


1.800 


0.410252 


20 


2.961675 


0.927128 


20 


2.15426 


0.556954 



Table 1 (L = 0) 



We presented here a fundamental (n = 0) overtone which is dominating 
in a signal. Some discussions of higher overtone behavior will be given in 
the following subsections. Note that in five space-time diamnsions we have 
a black hole-type solution only for a < 1.9, that is why the corresponding 
boxes in the table are empty. 

3.2 The neutral black hole, 1 — 1,2 and asymptotically high I 

The case of higher than zero multipoles is favorable for WKB methods. Here 
we observe good convergence of the QN values when increasing the WKB 
order. The I = 1, 2 fundamental QN modes are given in the tables 2 and 3 
respectively. 
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a [JJ = 5) 


i ~> 

Ke lo 


-1m oo 


a {V = 7) 


Ke lo 


-1m lo 


1/10 


0.720423 


0.255506 


1/10 


1.54573 


0.577608 


1/5 


0.723177 


0.252684 


1/5 


1.53194 


0.58701 


1/2 


0.739205 


0.236885 


1/2 


1.56277 


0.554551 


5 






5 


2.01379 


0.308316 


10 






10 


2.47824 


0.423737 


20 






20 


3.39094 


0.59452 


a (£> = 6) 


Re oo 


-Im oo 


a (D = 8) 


Re oo 


-Im lo 


1/10 


1.139007 


0.415034 


1/10 


1.93407 


0.750046 


1/5 


1.136193 


0.415706 


1/5 


1.90391 


0.773947 


1/2 


1.158635 


0.391339 


1/2 


1.94463 


0.734591 


5 


1.790103 


0.258382 


5 


2.44511 


0.455166 


10 


3.260415 


0.498426 


10 


2.6644 


0.463966 


20 


6.437139 


0.991374 


20 


3.2276 


0.581024 



Table 2 (L = 1) 



a {D = 5) 


Re lo 


-Im oo 


a {D = 7) 


Re oo 


-Im lo 


1/10 


1.07492 


0.248021 


1/10 


2.10314 


0.525815 


1/5 


1.08169 


0.243352 


1/5 


2.10106 


0.524304 


1/2 


1.10459 


0.228154 


1/2 


2.10473 


0.51385 


5 






5 


2.66996 


0.330935 


10 






10 


3.41504 


0.433789 


20 






20 


4.69346 


0.603341 


a (D = 6) 


Re lo 


-Im lo 


a (D = 8) 


Re oo 


-Im lo 


1/10 


1.60284 


0.39205 


1/10 


2.58782 


0.651764 


1/5 


1.607508 


0.387098 


1/5 


2.57296 


0.658318 


1/2 


1.62621 


0.370114 


1/2 


2.54771 


0.665021 


5 


2.619279 


0.269015 


5 


3.06761 


0.428616 


10 


4.81752 


0.505708 


10 


3.49737 


0.47698 


20 


9.52384 


1.00306 


20 


4.26009 


0.589931 



Table 3 (L = 2) 



For large multipole number / we can obtain the asymptotic formula using 
the first order WKB approach. For this to perform we first represent the 
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Re co 




a 



Figure 1: Real part of to as a function of the Gauss-Bonnet coupling constant 
a for I = 2, n = modes of six(diamond), seven(star) and eight(box) 
dimensional black holes 



-Im oj 




Figure 2: Imaginary part of w as a function of the Gauss-Bonnet cou- 
pling constant a for I = 2, n = modes of six(diamond), seven(star) and 
eight (box) dimensional black holes 
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value of r at which the effective potential attains the maximum in the form: 



W = r + n {l /I) + r 2 (l//) 2 + ... (13) 

Then we expand the effective potential and its second derivative in orders 
1 / 1 and insert this exapnsded values into the first order WKB formula. As a 
result for the most slow damping modes we obtain the asymptotic formula: 



-1 + Jl + 8rl- D aM 
L °Re = \— I 1 + A > UJlm = B ( 14 ) 

where A and B are constants which does not depend on I but depends 
only on black hole parameters and determined by cumbersome relations 
and therefore are not given here. 

To show that as a goes to zero the QN modes of GB black hole goes 
to their Schwarzschild values, let us look at the fit for some sample mode 
(D = 5, I = 2, n = 0), built on five values of a, a = 1/100, 1/200, ...1/500: 

uj = 1.0681 - 0.2528i+ 

(0.0422 + 0.0055i)a + (2.4021 + 3.8373i)a 2 + 0[a 3 ]. (15) 

We see that the fit is approching the corresponding Schwarzschild value 
1.0681 - 0.2529i. 

(Note that in [12] the tensor type gravitational modes are given which 
coincide with scalar field modes found in |4T)j). 

This result can be easily explained by the fact that for small a one has 
an asymptotic 

2M AaM ^ r 2l 

and therefore as a goes to zero the metric approaches the Schwarzschild 
metric. 

The dependence of the quasinormal modes on a can be seen from sam- 
ple figures 1 and 2 and all the above data. Namely we see that when a is 
vanishing, both real and imaginary parts approach its higher dimensional 
Schwarzschild values, while when a is becoming large, the quasinormal fre- 
quencies are proportional to a. Apparently it is possible to obtain some 
analytic formula for this regime of large a. Note that one can observe the 
convergence of the results obtained with the WKB formula when increasing 
the WKB order for both large and very small a. The worst convergence 
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occurs for intermediate values of a, yet this happens only for zero multipole 
modes when I = n. For I > n the convergence is good for the whole range 
of parameters. Note that the WKB method in the Will-Schutz form uses 
the series which is convergent only asymptotically, therefore we are unable 
estimate the error strictly within used method. Yet, in principle, it can be 
modified to control the error [IS]. All previous expiriance in usage of the 
Will-Schutz approach, based on comparison with accurate numerical results, 
says that if we observe the convergence, then it is the convergence to some 
true QN value, and the greater the WKB order the better the accuracy. Cer- 
tainly we do not have any garantee that we shall observe the convergence in 
any case we wish. 

Thus from (12) we see that for the GB black hole, similar to the higher 
dimensional Schwarzschild black hole behavior, for asymptotically large mul- 
tipole I, the real part is proportional to I while the imaginary part approaches 
some constant. The higher overtones within WKB method used here can 
be obtained only with I > n. For instance, from the following data for suf- 
ficiently large multipole (I = 5) we can reproduced first several overtones: 



4.73161 - 0.34419i n = 

4.70975 - 1.02805* n = 1 

4.64921 - 1.69859t n = 2 

4.49639 - 2.35764i n = 3 

4.13195 - 3.05609i n = 4 

In the above data a = 5, and D = 7. On this example data we see 
that the first few higher overtones have decreasing real part as the imagi- 
nary part increases. That feature is similar to the Schwarzschild black hole 
quasinormal behavior. To find sufficiently high overtones one should use the 
continued fraction technique. 



3.3 The charged black hole 

The dependence of the quasinormal spectrum on the charge of the GB black 
hole is quite different from the case of the Einstein black holes. In particular, 
from [11], [13], we know that the imaginary part of the quasinormal 
mode for Reissner-Nordstrem or dilaton black holes is growing as the charge 
is increasing until some its maximum value, then, usually at the charge 
equaled approximately to 0.8Q there is a sharp fall down. As can be seen 
from the Fig. 3, there is no such a maximum value of the damping rate, 
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Figure 3: Imaginary part of uj for D=6(diamond), 7 (star), 8(box) dimen- 
sional black holes on example of L = 2, n = 0, a = 5. 



but instead, the monotonic decreasing exists. On the Fig.4 one can see that 
the real oscillation frequency is just monotonically increases as the charge 
of the black hole is increased. 

At high I and non extremal value of the black hole charge, in a similar 
fashion to the case of the neutral black hole, we obtain the asymptotic 
formula: 



-1 + y/l + 8rl~ D aM - gggg 



1/2 



^ e =|^ U 2q 5z"tt^| l + A >, U Im = B>, 

(16) 

where the constants A' and B' reduces to A and B when Q = 0. The explicit 
form of these constants and of the coefficient ro, ri..etc. is available form 
the author upon request. 
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Figure 4: Real part of w for D=6(diamond), 7 (star), 8(box) dimensional 
black holes on example of L = 2, n = 0, a = 5 
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4 Discussions 



We have investigated quasinormal modes of a black hole in Gauss-Bonnet 
gravity. It is shown that quasinormal behavior is totally different from that 
of black hole in Einstein gravity. In particular, the quasinoraml spectrum 
essetailly depends on the Gauss-Bonnet coupling parameter: at large a the 
QN modes are proportional to a, while when a goes to zero, the QN modes 
approach their Schwarzschild values. Due to the Gauss-Bonnet coupling the 
quasinormal modes have greater oscillation frequency and greater damping 
rate at large a, but at moderate a, the damping rate as a function of a, first, 
is decreasing until some minimum value and then begin to grow by the low 
~ a. As can be seen from comparison of the scalar filed potential considered 
here with that found in [T3|, for tensor-type gravitational perturbations, the 
scalar field potential and the tensor type gravitational potential coincide. 
The same situation takes place for higher dimensional black holes in Ein- 
stein gravity [Ha]- Thus the quasinormal spectra are the same for adequate 
multipoles. If to disentangle equations for other types of gravitational per- 
turbations of GB black hole is a solvable problem (see for instance references 
in |46| . |47j ) . then the question of stability and gravitational QN radiation 
would be completed. 

The problem that was beyond the present study is the finding of the 
asymptotically high overtones of the QN spectrum. After all, the quasinor- 
mal modes for the GB black hole in Anti-de Sitter space-time may have the 
ADC/CFT holographic interpretation and are also interesting to study. We 
hope that further investigation will clarify these points. 
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